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Abstract 

In the deconfined phase of quark-gluon plasma, it seems that most of the quarks, anti- 
quarks and gluons should be effectively free in the absence of the linear confining po- 
tential. However, the remaining Coulomb-type potential between quarks in the plasma 
could still be sufficiently strong that certain bound states, notably of heavy quarks such 
as J/ip are stable even in the deconfined plasma up to a certain temperature. Baryons 
can also exist in the deconfined phase provided that the density is sufficiently large. We 
study three kinds of exotic multi-quark bound states in the deconfined phase of quark- 
gluon plasma from gravity dual models in addition to the normal baryon. They are 
/c-baryon, [N + /c)-baryon and a bound state of j mesons which we call "j'-mesonance" . 
Binding energies and screening lengths of these exotic states are studied and are found 
to have similar properties to those of mesons and baryons at the leading order. Phase 
diagram for the exotic nuclear phases is subsequently studied in the Sakai-Sugimoto 
model. Even though the exotics are less stable than normal baryons, in the region of 
high chemical potential and low temperature, they are more stable thermodynamically 
than the vacuum and chiral-symmetric quark-gluon plasma phases (xS-QGP). 
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1 Introduction 



The discovery of AdS/CFT correspondence [TJ [2] provides a new tool for studying the 
strongly coupled gauge theories. Although the original setup and most of the systems that 
string theorists have been investigating so far are highly supersymmetric and conformal, a 
lot of progress has been made in constructing more realistic models. Now we have exam- 
ples of QCD-like gauge theory with known gravity dual that share most of the qualitative 
features of QCD. These holographic models allow us to perform analytic calculations in the 
regimes which are too difficult to implement for the real QCD even for lattice calculations. 
The properties of quark-gluon plasma from Relativistic Heavy Ion Collisions and QCD at 
finite baryon density are two examples of such regimes. 

The gravity dual of baryons can be described via baryon vertex [SI H] , a D-brane wrapping 
higher dimensional sphere in 10-dimensional curved background with N strings attached to 
it and ending at the boundary. These strings are required to cancel an N charge in the 
world- volume of the wrapped brane due to the presence of RR flux in the background. The 
endpoint of fundamental string that ends on D-brane is electrically charged with respect 
to world-volume U(l) gauge field. Its charge is +1 or —1 depending on the orientation 
of the string and D-brane. Moreover, strings stretching from the baryon vertex to the 
boundary of AdS or the corresponding background spacetime (e.g. in Sakai-Sugimoto model) 
behave as fermions, giving antisymmetricity to the baryon vertex. This fact allows us to 
construct an SU(N) gauge-invariant combination of N quarks as required by the group 
theory. Baryon configurations were investigated further in [5]- [7]. The authors in [5] extended 
the consideration in confining background where it was found that the binding energy is linear 
in N and in the size of the baryon on the boundary. And furthermore, they found that in 
Nsusy = 4 theory there are stable configurations for baryons which are made of k quarks, 
or "k-baryon", if 5iV/8 < k < N. Such configurations can be realized by considering the 
usual baryon vertex with k strings stretched up to the boundary and the rest N — k strings 
stretched down to the horizon. These baryons are not colour singlet and transform as k \^i k y 
representation under SU (N) gauge group, for example the case k — N — 1 gives rise to a 
baryonic configuration in the anti-fundamental representation. In a confining theory we do 
not expect to find such a bound state. It was proposed in [9j [10] that the k < N bound 
states can exist in a deconfined phase. 

In general, we could imagine that there would be more exotic baryon states in the decon- 
fined phase where bound states of quarks need not be the colour singlet. Some attempts have 
been made in constructing holographic description of exotic multi-quarks bound states [9]- 
[T2] . The author in [12] considered exotic quark configurations formed by combining two or 
more baryon vertices together. However, it might be possible to construct an exotic baryon 
from a single baryon vertex which should be more energetically preferable. One useful ob- 
servation is that there are infinite combinations of string charges that can cancel the charge 
from the background RR flux. Hence, the total number of strings attached to the baryon 
vertex need not to be equal to N. For example, if the orientation of D-branes is fixed in such 
a way that there is +N units of charge on its world-volume, we can attach N + k strings, 
each with —1 charge and k strings with +1 charge to make the total charge vanishes. As 
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long as the conservation of charge is concerned, k could be any integer. In this case, we can 
construct a k > N baryon. Such baryon could be the lightest bound state in some irreducible 
representation of the underlying gauge theory thus it may be stable and can be observed in 
the deconfined phase. We would like to investigate this possibility further in this paper. 

It is also interesting to study exotic baryons in more realistic model such as Sakai- 
Sugimoto model [131 E]. This model is based on Witten's model [15] which uses the D4- 
brane wrapping a Scherk-Schwarz circle and adds a stack of Nf probe D8-branes and a stack 
of Nf probe anti-D8-branes transverse to the circle. This model contains massless chiral 
fermions and the flavour symmetry. The most striking feature of this model is that it in- 
troduces geometrical mechanism for spontaneous chiral symmetry breaking. Using the fact 
that the circle vanishes at a finite radial coordinates in the near horizon limit, D8-branes 
and anti-D8-branes are connected in a U-shaped configuration. At low temperature the 
model describes a confining gauge theory with broken chiral symmetry. Above a deconfine- 
ment temperature, gluons become effectively free. However, both the connected U-shape 
D8-branes configuration and the separated parallel brane-anti-brane configuration are pos- 
sible in the intermediate temperature. The chiral symmetry is still broken even though the 
gluons are already deconfined. At higher temperature the chiral symmetry is restored, which 
is illustrated geometrically by the separation of the D8-branes and anti-D8-branes [IB] . This 
corresponds to the branes being in parallel configuration. 

The model also has an interesting phase structure. Finite baryon density in the Sakai- 
Sugimoto model has been studied in [T71 [18] and extended to the full parameter space in [T9] 
where baryon matter is represented by D4-branes in the D8-brane (nuclear matters) or by 
strings stretched from the D8-brane down to the horizon (quark matters). It was shown that 
the former configuration is always preferred to the latter and quark matters are unstable 
to density fluctuations. In the deconfined phase there are three regions corresponding to 
the vacuum, quark-gluon plasma, and nuclear matter, with a first-order and a second-order 
phase transition separating these three phases. The author in [TH] found that for a large 
baryon number density, and at low temperatures, the dominant phase has broken chiral 
symmetry in agreement with QCD at high density. It is interesting to see how exotic baryon 
states fit into the phase structure. 

This paper is organized as the following. In section 2, we discuss some classes of exotic 
baryon configurations and investigate their static configurations in section 3. Binding energy 
and screening length of the configurations are calculated in section 4. The dependence on 
free quark mass of exotic baryon configuration is discussed in section 5. The phase diagram 
of Sakai-Sukimoto model with exotic baryons is investigated in section 6. We discuss our 
results in section 7 and conclude in section 8. 

2 Some classes of multi-quark states 

In the deconfined phase of QGP, coloured states of a number of quarks and antiquarks can 
exist in the medium as long as it is energetically more favoured than the free quarks and 
antiquarks or other mesonic states. We will call these multi-quark states as "baryons" in 
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this article. In the confined phase, the only allowed baryons are those with colour singlet 
combinations such as nucleons and pentaquarks. For the deconfined phase, baryons can have 
colour and thus can have more varieties than the situation in the confined phase. 

In general, a D(8 — p)-brane wrapping the subspace S$- p of the background spacetime 
sources the gauge field An\ on its world volume. This gauge field will couple with the 
antisymmetric (8 — p)-form field strength GVs-p) and induce the charge upon the wrapping 
D(8 — p)-brane. If the background is generated by a stack of N Dp-branes, then the charge 
being induced upon the wrapping D(8 — p)-brane will be exactly N. This charge needs to 
be cancelled by external charges brought about by strings. Each of the strings stretching 
out from the wrapping brane to the spacetime boundary or probe branes carries —1 unit of 
charge. Therefore it is required that the total number of "quark" strings stretching out from 
the wrapping brane must be N. The configuration of wrapping D(8 — p)-brane with totally 
N strings stretching out is called a baryon vertex [31 0] . 

For the confined phase, since quarks cannot exist as free-quark strings with one end 
falling behind the horizon, therefore they can only start from the baryon vertex and go to 
the probe branes. On the other hand, in the deconfined phase, a radial string configuration 
lying along the radial coordinate is also a classical solution of the Nambu-Goto action [2U] 
and it represents the free (anti)quark state in the QGP medium. A string can either start 
from the baryon vertex and go radially to the horizon of the background spacetime or it can 
come from the horizon and end at the baryon vertex. We will call this string configuration 
which is allowed in the deconfined phase as the "radial string" . 

In the deconfined phase of QGP, it is possible to have strings hanging from the space- 
time boundary down to the baryon vertex and another k T strings stretching radially from the 
baryon vertex down to the horizon. The total number fch + k T = N is the charge conservation 
constraint on the configuration. This configuration is known as "A>baryon" [8]. 

Another possible configuration is when there are N quark-strings and k antiquark-strings 
hanging down to the vertex from the probe branes. To conserve the charge, there are 
additional k quark-strings hanging from the vertex down to the horizon. We will call this 
configuration "(A r + fc)-baryon" (e.g. pentaquark could be represented by one of this kind). 

An even more interesting configuration allowed in the deconfined phase is when there are 
j pairs of quark and antiquark strings hanging from the probe branes down to the vertex. 
Again, to conserve charges, we need N radial strings stretching from the vertex down to the 
horizon. This configuration obviously can decay into j mesons when it is less energetically 
favoured. Therefore we will call this state, a "j-mesonance" , representing a binding state 
of j mesons in the QGP. 

In summary, the charge conservation constraint for each case can be expressed as the 
following. 

For /c-baryon, 

k h + k r = N; k h = k. (1) 

For (N + /c)-baryon, 

k h -k r = N; k h = N + k. (2) 
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Figure 1: The gravity dual configurations of the hypothetical exotic states (a) /c-baryon with 
the number of hanging strings k^ = k < N and the number of radial strings k T = N — k. (b) 
(N + fc)-baryon with k^ = N + k and k T = k. (c) j-mesonance with k^ = 2j and k T = N. 



For j-mesonance, 

k h = 2j; k r = N. (3) 

Note that k^ is the number of strings hanging from the boundary down to the baryon vertex 
and k T is the number of strings hanging from the vertex down to the horizon. The value of 
k and j can be as large as N x Nf. However, in this article, we will take this number to 
be large and ignore the upper bound on k and j. Each configuration of exotic baryons is 
illustrated in Fig. [TJ 

3 Force conditions 

In this section, we will consider the force condition for each exotic configuration of the 
quarks and antiquarks in a deconfined phase. The calculation will be performed in the 
gravity background similar to those of Sakai and Sugimoto's [13]. Even though the chiral 
symmetry restoration can be addressed within this model, we will not consider the aspect in 
this section but rather focus our attention on the high temperature phase where quarks and 
antiquarks are effectively free in the absence of the linear confining potential. The positions 
of D8/D8 will be taken to be large and we will approximate it to be infinity in this section 
as well as in the discussion of binding energy and screening length in section HI Analysis in 
this heavy-quark limit provides us with valuable physical understanding of certain essential 
features of the exotic states. Generalized results for a near-horizon background metric of 
the Dp-branes solution and its dependence on positions of the probe branes will be given in 
section 

Even in the deconfined phase, quarks and antiquarks feel effective (screened) potential 
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from other constituents. Therefore, a number of population of them will exist in various 
forms of bound states, some of which are exotic in the sense that they cannot be formed in 
the confined phase at low temperature. 

Start with the following background metric 



\RdaJ ' \ u J V /(i 



F (4) = ^e 4 , e * = g s (-^-) , R 3 m = 7ig s Nll 

where f{u) = 1 — u^/ii 3 , ut = I611 2 R^T 2 / 9 . Note that the compactified £4 coordinate (x 4 
transverse to the probe D8 branes), with arbitrary periodicity 2ttR, never shrinks to zero. 
The volume of the unit four-sphere Q4 is denoted by V4 and the corresponding volume 4-form 
by e 4 . F(4) is the 4-form field strength, l s is the string length and g s is the string coupling. 
The dilaton in this background has w-dependence and its value changes along the radial 
direction u. This is a crucial difference in comparison to the AdS-Schwarzschild metric case 
where dilaton contribution is constant. 

The action of the baryon configuration is given by 

S = S m + k h S F1 + k T S F1 , (4) 

where Sda represents the action of the D4-brane. Spi is the action of a stretched string from 
the boundary down to the baryon vertex and Sfi is the action of a radial string hanging 
from the baryon vertex down to the horizon. Recall that Sda can be obtained from the 
Dirac-Born-Infeld action^- After some calculations, we obtain 



where r is the total time over which we evaluate the action and u c is the position where the 
D4-brane vertex is located. 

The variation of the action with respect to u gives the volume term and the surface term. 
The volume term leads to the usual Euler-Lagrange equation for the classical configuration of 
strings. As an approximation, we assume the baryon vertex to be a point (not being distorted 
by the connecting strings) located at a fixed value of u = u c as in Ref. [5|. Under this 
assumption, the surface terms provide additional zero-force condition on the configuration, 

jG (x) - KB + k r = (6) 



S DB i = J dx°dt P T p ; T p = (e^(27r)Pa' (p+1)/2 ) ' y/-det(g) 
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where 



1 i / 
1 + tt u T it 



G (x) = 2 , s = — < 1, and B= c . (7) 

vi 3 ^ «« v/ <2 + /K)( ^ )3 

Notice that these conditions occur at the location of the vertex at u — u c , at which there 
exists the balance between the pull-up force (toward the direction of increasing u) due to 
the tension of hanging strings and the pull-down force due to the "weig ht"d of D4-brane plus 
the tension of radial strings. 
Since B < 1, we obtain 

N 

k h >jG (x) + k r , (8) 

which expresses the lower bound of the number of hanging strings. In other words, the 
number of hanging strings cannot be less than this critical value, otherwise the no-force 
condition is not satisfied. The equality of ([8]) is held only when all hanging strings are 
stretched straight, otherwise we require more hanging strings to balance the pull-down force. 
Let us now consider each class of the multi-quark states. 

In the case of /c-baryon, plugging the condition (TjQ) into (JH1), we obtain 

K = k> ^(G ( x ) + 3). (9) 
o 

Apart from the lower bound, we also have the upper bound, k < N, therefore Gq(x) cannot 
be larger than 3, resulting in 

x< 0.922. (10) 

Notice that this restriction on x is a result from the conditions of the force balance and 
conservation of string charges. This shows that there is an upper-bound on the temperature, 
over which the horizon is too near to the point vertex that the pull-down force always 
overcomes the pull-up one. 

In the case of (A" + /c)-baryon, in the same way as the preceding case, plugging the 
condition of charge conservation (T5]) into ([8]), we have the following condition, 

N 

k h = N + k> —G {x) + k. 
3 

Unlike the case of fc-baryon, the upper-bound of the number of hanging strings does not 
exist. However, we still obtain the same condition Gq(x) < 3, hence x < 0.922. 
Finally, in the case of j-mesonance, similarly, Eqn. ([3]) results in 

J>^(G (x) + 3). (11) 
o 



2 This is not exactly the weight in the usual sense since the direct gravitational force on Dbrane is already 
balanced by the force from the RR-flux, but it is the force originated from minimization of self-energy due 
to the brane tension caused by the background metric and the gauge interaction. This is very similar to 
the self-energy of a spring under gravity where the spring potential energy changes with the tidal force from 
gravity in the background. The DBI action of the D4~ u c \J f{u c ) which is positive for u c > ut and becomes 
zero (minimum) at u c = ut and thus it represents the "weight" on D4 towards the horizon. 
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The lower-bound of the value of j is 2N/3 at zero temperature (x = 0) and it will be larger 
as the temperature grows. Nevertheless, the upper-bound of the limit on j does not exist. 

Finally, we would like to comment on the limits on the value of k, k,j when the temper- 
ature is zero. In terms of n = 7 — p (of the spacetime background generated by Dp-branes), 
the condition (jSJ) becomes 

h > — + K (12) 

n 

which leads to 

1 L > 2+1 (i3) 

and no conditions on k. This critical numbers are 5/8, 2/3 for n = 4, 3 (the AdS-Schwarzschild 
and Sakai-Sugimoto model) respectively. It is an interesting coincidence that the critical 
numbers are the same for both /c-baryon and j-mesonance. Even though it appears from 
Eqn. (j!2p that there should also be a constraint on the (N + k) configuration, it turns out 
that there is none. 



4 Binding energy and the screening length 

In this section we will calculate the binding energies of the /c-baryon, (N + /c)-baryon, and 
j-mesonance in the deconfined phase. These binding energies are taken to be the differences 
between the total energies of each configuration and the corresponding energies of the free 
strings configuration which represents the free quarks and/or antiquarks state. The number 
of free strings in the free quarks state is determined solely by the total number of strings 
hanging from the boundary, kh- 

The total energy of each configuration is given by E = S/r of the corresponding action 

5 for each configuration. The binding energy for each hanging string is consequently, 




\ 3 

u \ „ , s 1 



( /(«) - ^ J du. (14) 

Due to the no- force condition in the surface term, we impose Eqn. (jSJ) and Eqn. (JTj), or 

where the tension of each hanging string at u c is constrained by 

/V k 

B = B(k h , kr, x) = — G (x) + f . (16) 
3fc h k h 

Since the Lagrangian C does not depend on a explicitly, the conserved Hamiltonian can 
be defined to be 

dC 

H = £ — u' — — = const, (17) 
du' 
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leading to 



/K)(t^) 3 



/H(^) 3 



V u c 2 + /K)(^i) 3 v /m ' 2 + /(m)( ^ )3 ' 

Then substituting Eqn. (|T5|) into this equation, we obtain 



;is) 



u 



/(«) 



/( Uc )(^-)3(l_ jB 2) — Vi?D4 

This gives the size (radius) of the baryon as seen on the gauge theory side, 

ul 12 Ji y V (V 3 - * 3 )(y 3 - ar* - (1 - x 3 )(l - B^)) " 

Note that u c « -^J 1 at the leading order. 

Using Eqn. (119j) and let y = u/u c , the regulated binding energy now becomes 



(19) 



(20) 



E 



u 



Fl 



(y 3 -x 3 ) - (1 -a; 3 )(l -S 2 ) 



(21) 



Hence, we obtain the total energy of the configurations as 



E 



Nu T ( Vl 



2tt 

AT 2 



3x 



N ) x \N 



x 



X 



(22) 
(23) 



where £ represents the terms within the brace of (I2TT) . 

To obtain the relations between the total energy of the configurations E(x) and L(x), 
we eliminate the parameter x = ut/u c . By numerical calculations, the results are shown 
in Fig. I2ll3l The binding energy of iV-baryon is the deepest, suggesting that it is the most 
tightly bound state. For (N + fc)-baryon, increasing k makes the binding energy smaller and 
the bound state is less tightly bound. The case of j-mesonance is quite similar. Generically, 
a j-mesonance has shallower binding potential than the total energy of j mesons. However, 
as j grows, the difference gets smaller and smaller. 

The screening radius or screening length of exotic multi-quark state is defined to be the 
value of radius L* at which the binding energy becomes zero from negative values at smaller 
distances. This screening radius is therefore one-half of the usual definition of screening 
length in the discussion of mesonic state where it is defined as the zero-potential distance 
between quark and antiquark. 

Numerical results suggest that the screening length of baryons and mesonance decrease 
as the temperature increases, i.e. L* ~ l/T for a fixed value of k, k, j as is shown in Fig. QMl 
This is the generic form for the screening length in both the AdS-Schwarzschild and Sakai- 
Sugimoto models because it is the quantity which does not depend on the 't Hooft coupling 
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Figure 3: Comparison of the potential per N between j-mesonance and j mesons for ji/N — 
0.8, ji/N = 3 at temperature T = 0.25. 
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Figure 4: Screening length with respect to k for the temperatures in 0.15 — 0.35 range. 
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at the leading order [21]. It is also an increasing function of k and j. Interestingly, (N + k)- 
baryon has the opposite tendency with the screening length decreases as k grows. On the 
other hand, the screening length of j-mesonance has a saturation value L*_ mesonance 
as j — > oo. 



r * 

meson 



5 Dependence on the free quark mass 

In this section, we will study dependence of the binding potential on the position of the 
probe branes. This is useful when position of the probe branes are at finite distance from 
the black hole horizon and the corresponding quarks have finite mass. For example, the 
probe branes are D8 and D8 flavour branes in the Sakai-Sugimoto model. 

The calculation of binding energy as a function of the radius L of the multi-quark states in 
the previous sections can be generalized to the case where the background metric is generated 
by a stack of Dp-branes as the following. Start with the energy of a hanging fundamental 
string with n = 7 — p, 



E 



u 



Fl 



7T- dy 



yti _ x r 



(y n - x n ) - (1 - x n )(l - A(n) ) 



x) 



(24) 
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and the radius, 



R n/2 



dy 



[1 - x n )(l - A(n) 2 ) 



Uc (n-2)/2 J i y (yn _ x n^ y n _ x n - (1 - X n )(l - A(n) 2 )) ' 

The total regulated binding energy of the configuration then becomes 

Nuh I y/l — x n , ( k\\ £ , f k r \ 1 — x 



(25) 



E, 



tot 



2n 



nx 



N x 



N I x 



(26) 



where 



£ 



dy 



yU 



x' 



(y n - x n ) - (1 - x n )(l - A(n) ) 



and 



A(n) = 

The parameter x is again given by 

u T {n) 



< 2 + /K)(7^) n 



N 
nkh 



1 + ^x n 

y/1 - X n 



(27) 



(28) 



.r 



= 3, 4) 
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7T 2 i? 3 T 2 , 7Ti? 2 T. 



(29) 



Note that the case n = 3 and n = 4 corresponds to the case of Sakai-Sugimoto and AdS- 
Schwarzschild gravity dual model respectively. 

Introduction of quark masses into the configuration can be done by terminating hanging 
strings at certain radial distance u max < oo. The universal behaviour of heavy-quark poten- 
tial comes from the limit u max — > oo. We can split the total binding potential of the string 
into two parts. The first part is the binding potential in the u max — > oo limit and the second 
part is the mass dependent potential. Therefore, the mass dependence part of the binding 
potential, EFi{u max ) (m = u max /27i), can be expressed as 



i?Fi (finite mass) 

U"max) 



(U r 



""max 
1*00 



Uc_ 

Urn 



00) + -E'fi( Umax ) 1 
dy 



yn _ x n 



Umax I UC 



; (1- A{n) 2 ) fu n c -u n 



Air(n - 1) 



(y n - x n ) - (1 - x n )(l - A(n) ) 

n y 

' ok 



- 1 



u. 



l-2 n] 
"max J 



(30) 
(31) 

(32) 



Eliminate u c by using 

R n/2 



L 



u (n-2)/2 



Umax I Uc 



dy 



x n )(l - A{nf) 



(y n - x n )(y n - x n - (1 - x n )(l - A(n) 2 ))' 



(33) 
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The result involves complicated functions of A which can be cast in the following form, 

E FX {u max ) ~ -ul^ x (R n2 l^h{A) + u n T f 2 {A)), (34) 

where fi t i(A) are some functions of A. 

Interestingly, the mass dependence of multiquark potentials has similar form as the mass 
dependence of mesonic state ~ m 1_n in Ref. [20]. This is natural due to the fact that most of 
the mass of constituent quarks come from the tail part of strings which extend to the large- 
u region. The mass dependence of the binding potential at the leading order is therefore 
determined only by the contribution of the hanging strings from the large-u region. As long 
as the background spacetime of the gravity dual is asymptotically similar to the background 
considered here in the large- u limit, we would expect the same mass dependence as the form 
we obtained in this section. 

6 Phase diagram 

A natural question to ask is whether we have a phase where the exotic multiquark states 
are preferred over the normal nuclear matter (namely the gas of iV-baryons), vacuum, and 
the chiral-symmetric quark-gluon plasma phase. To consider a realistic model where these 
three phases are distinct, we focus our consideration on the Sakai-Sugimoto model with 
n — 3. To calculate the phase diagram involving exotic states, it is necessary to consider 
the contribution from D8 and D8-branes in the Sakai-Sugimoto model in addition to the 
contributions from strings and D4-branes. We will assume that the characteristic distance 
between D8 and D8 in x 4 direction is Lq. The relevant scales of the model therefore depend 
on both u T and L . 

When there is no radial string pulling the vertex down towards the horizon, it was 
demonstrated in Ref. [7] by numerical method that the vertex will be pulled all the way 
up to the position of the flavour branes if the temperature is not very high. Addition of 
radial strings to the vertex would pull the vertex and the flavour branes towards the horizon. 
As temperature rises, the radial strings pull the vertex down with stronger force since they 
are closer to the horizon. It is possible that the vertex then starts to separate from the 
flavour branes and we might need to consider the configuration where the vertex and flavour 
branes are separated. However, we can see that the difference between the two configurations 
should be relatively small (namely, only the force conditions will be slightly different) and 
we should be able to approximate the situation by considering the configuration where the 
vertex is not separated from the flavour branes. It is also shown in the Appendix that this 
configuration satisfies the force condition and thus is allowed. Therefore, it will be assumed 
that the vertex is always in the flavour branes for the discussion in this section. Moreover, 
the vertex will be treated as a static configuration and any distortion caused by the strings 
attached to it will be ignored. 

The calculations presented in this section are adapted from Ref. [TU] except that we add 
radial strings hanging from the vertex down to the horizon for the consideration of exotic 
nuclear phase. We also use position of the D4, u c , instead of uq (where ^(uq) — > oo) in our 
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calculation concerning the exotics. This approach allows us to deal with the contribution 
from radial strings more conveniently. As is shown in Fig. the vacuum phase with broken 
chiral symmetry corresponds to the configuration where D8 and D8 are connected into a 
curve in the x<± — u projection. The chiral-symmetric phase of quark-gluon plasma (%S- 
QGP) corresponds to the configuration with the parallel D8 and D8 stretching from the 
spacetime boundary down to the horizon. Finally, the nuclear (including exotics) phase 
corresponds to the configuration where the D4 vertex is located at the D8-D8 curve, pulling 
it down towards the horizon by its "weight" in the background. Each vertex has radial 
strings attached to it, pulling it further towards the horizon. When there is no radial strings 
attached, the nuclear phase is of normal iV-baryons. The chiral symmetry is also broken in 
this phase. 

Under the above assumptions, the contribution from the strings hanging down from the 
spacetime boundary to the vertex is negligible. The only contribution of strings is from 
the radial strings hanging down from the vertex to the horizon. The total action of the 
configuration is given by 

S = Sr>8 + Sd4 + Spi- (35) 
Generically, the DBI action of the D8-branes is given by 

Sds = -/is J d 9 Xe^Ti^-det(g M N + 2na'F MN ) (36) 

where the field strength of the flavour group U(Nf) is 

T = dA + %A A A. (37) 

The flavour branes provide "global" quantum numbers such as baryon number to the string 
and subsequently to the strings-brane configuration dual to baryon in the gauge theory side. 
The diagonal part of the representation matrix of U(Nf) is the U(l) subgroup which induces 
baryon number to the end of string attached to the flavour branes. Redefine the U(l) part 
so that 

A = A su{Nf) + —=A (38) 
V 2iV / 

with A represents the U(l) piece of the gauge field. The DBI action of the D8-brane coupled 
to the diagonal gauge field is then given by 



S D8 = M J du mV/(«)K(«)) 2 + «" 3 (1 - WW) ( 39 ) 
where the constant scales linearly with Nf as 

9s 
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and the rescaled U(l) diagonal field, 



2na'A 



(41) 



The action does not depend on do(u) explicitly and therefore a constant of motion can be 
defined as 



d 



ua {u) 



V/(«)K(«)) 2 + «- 3 (i-K(«))- 



(42) 



We will see below that the constant d can be interpreted as the baryon number density 
sourced by the D4-branes once we introduce the Chern-Simon action of the gauge field. 
Note that d plays the role of the electric displacement field [19]. In the confined phase, the 
only possible source for d is the D4-brane wrapped on S A in the D8-branes. In the deconfined 
phase, either D4-brane or strings which stretch from the D8-brane down to the horizon can 
serve as the sources for d. Here, in the study of exotic baryons, we consider the case where 
both D4-brane and strings are present as the sources. This possibility was not investigated 
in p2]. 

Similarly, the constant of motion with respect to x±(u) leads to 



1 



f(u)(u 8 + u 3 d 2 ' 



u 3 f(u) lf(u )(u 8 + u 3 d 2 ) 



(43) 



where u is the position when x' 4 (u ) = oo. 

Instead of using w as the reference position, the radial position of the D4 on the D8- 
branes, u c , can be used to calculate x' 4 (u), 



where 
F 



f(u c )^u 8 c + u 3 c d 2 
\/ f{u c )(x' A {u c )) 2 + w c 



u 3 f(u) 



f(u)(u 8 + u 3 d 2 ' 
F 2 



-l 



(44) 



(45) 



V u3 cf( U c) 



%ul + d 2 ) 



d 2 



i + |(^) 3 + 3^V7( 



u. 



f(Uc 



,(46) 



The number of radial strings n s represents the number of strings hanging down from D4- 
branes to the horizon in unit of 1/N. For k, (N + fc)-baryon and j-mesonance, the values of 
n s are 1 — k/N,k/N,l respectively. Calculation of x' A (u c ) is performed by minimizing the 
action with respect to the variation of u c (see Appendix). For a fixed Lq, increasing the 
number of strings n s results in D4-D8 configuration being pulled down more towards the 
horizon. 
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The U (Nf) gauge field A also generates Chern-Simon term, 



Sea = / <»*(A). (47) 

24tt z J M 4 xR 

For A = A^dx^ 1 + A u du, the 5-form field strength is given by 

uM) = Tr(AF 2 -^A 3 F+^A 5 y (48) 

Only the first term contains non-vanishing contribution from the U(l) part which would 
be identified with the number density of baryon. We will assume a uniform distribution 
of the gas of D4-branes in M 3 at u = u c in the radial direction. This leads to the relation 
between the number density of D4-branes, 714, and baryon number density d [T9] . 



n * = ^r d - (49) 

Phase transition for a system where the number of particles varies is most conveniently 
described by the grand canonical ensemble. The grand canonical potential of each phase can 
be defined using the corresponding action of the D8-branes as 

fl(/x) = ^S D8 [x 4 (u),do(u)} d . (50) 

The baryon chemical potential is given by the U(l) diagonal field at the boundary, 

/i = a (00), (51) 
from which the baryon number density is determined, 

d = (52) 

This justifies the association of grand canonical potential with the D8 action. When addi- 
tional sources of the baryon number are introduced, the free energy, JF E , from the sources 
will also contribute to the baryon chemical potential, 

d 1 / ~ \ dJ~ 

H = QjTf (J=>D8[x4(u),d(u)] d + S source (d,u c )J = -Qj-, (53) 

where the Legendre-transformed action Sds is given by 

POO 

Sds = S D8 +N d(u)d' du, (54) 



M f°° du uV/(u)«(u)) 2 + n" 3 J 1 + ^. (55) 
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In our case, the additional sources are D4 and radial strings. These relations can also be 
applied to the vacuum phase (with broken chiral symmetry) where u c is replaced with w . 

Setting Lq = 2 J^° =UQU x' A {u)du = 1, the expressions for the grand canonical potential 
and the chemical potential for each phase are given by 

vacuum phase, d = 0: 



xS-QGP phase, 



Xa[U\ 



0: 



du 



u 



5/2 



Uq 



/(«) - 5/(«o) 



(56) 



Hqgp 



du—z 

UT VM 5 + d? 

00 d 
du—z 
UT VM 5 + d? 



(57) 
(58) 



nuclear (including exotics) phase : 



du 



Uc 
OO 



f(u)(u 8 + uW)l ^FTd? 



-1/2 



U ' 



(59) 



du 



-1/2 



/ t— = + o^V 7 7K) + w«(m c - Mr)- (60) 



At a fixed temperature T and chemical potential /x, a first order phase transition line between 
phase 1 and 2 is obtained when Qi = f2 2 ,/ii = /i2 = A 4 - Transitions between vacuum <-> %S- 
QGP and xS-QGP <-> nuclear phases are of this kind. On the other hand, phase transition 
between nuclear «-> vacuum is second order in nature, at least for this case when there is no 
interaction between each D4. The second order phase transition line occurs when 



dd 



d 2 r E 
dd 2 



(61) 



has discontinuity at d — 0. 

In the Sakai-Sugimoto model, there is a phase transition temperature above which glu- 
ons become deconfined. However, it does not necessarily imply that everything including 
quark and antiquark is totally free and chiral symmetry is completely restored above this 
temperature. When the baryon chemical potential is sufficiently high, baryons can exist 
even when the temperature is higher than the deconfinement temperature [19]. Only when 
the temperature increases even further that everything will be completely dissolved and the 
chiral symmetry is also restored. We also see this behavior in the phase diagram in Figure [8] 
where we ignore the confined region at low temperature and present only the deconfined 
part of the phase diagram. 

The phase diagram of vacuum with broken chiral symmetry, xS-QGP and phase of nu- 
clear including exotic multi-quark states is shown in Figure [H The phase diagram involving 
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(a) (b) (c) 



Figure 7: Configurations of xS-QGP (a), vacuum (b) and exotic nuclear phase (c) in x A — u 
projection. 



0.175 - 




Figure 8: The phase diagram of exotic nuclear matters above the deconfinement temperature. 
Nuclear phase including exotics is shown as the region on the lower right corner where 
it is divided into 3 parts for representative purpose. A, B, C represents the region where 
exotic baryon phase with n s = (iV-baryon), 0.1, 0.3 is preferred over vacuum and xS-QGP 
respectively. 
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vacuum and xS-QGP phases was first obtained in Ref. [IS] and the full phase diagram 
without the exotics was obtained in Ref. [TP]. Since the strings pull down the D4-D8 con- 
figuration towards the horizon, the configuration with n s > is less stable than the normal 
iV-baryon (n s = 0). This is shown in Fig. [8] where the region of n s > nuclear phase (B, C) 
is smaller than the region of iV-baryon phase (A). They are actually less stable than the 
iV-baryon since the grand canonical potential Q ng>0 (T, fi) > Q ng=0 (T, /j,) for 0.5 > n s > 0. 
Above n s > 0.3, the exotic phase becomes unstable to density fluctuations (|^ < 0) at high 
temperatures in certain range of d but still remains stable in a region of parameter space. 
Numerical studies reveal that for approximately n s > 0.5, the multiquark states become un- 
stable thermodynamically with respect to density fluctuations for most of the temperatures. 

Addition of radial strings introduces extra source of the baryonic chemical potential. 
We can see from Fig. [S] that the value of fi on set for the exotic nuclear phase increases with 
the value of n s . Nevertheless, once emerged (i.e. /i > /i onse t), the exotic phases are more 
stable than the vacuum at any temperature, but less stable than xS-QGP at sufficiently high 
temperatures above which chiral symmetry is restored. 



7 Discussions 

It is desirable to compare the binding energy of each multi-quark state in order to discuss the 
stability of each configuration as well as their relative abundances in the deconfined phase. 
At a fixed temperature T, we can compare numerically the binding energies E as functions 
of the size L of the configuration as is shown in Fig. [2131 For /c-baryon and (N + /c)-baryon, 
we compare the energy with iV-baryon. For j-mesonance, we compare the energy with the 
energy of j mesons. 

From Fig. [21 iV-baryon is more energetically favoured than /c-baryon and (N + /c)-baryon 
for any value of k, k. Since there are less hanging strings from the spacetime boundary and 
more radial strings pulled down into the horizon in the case of /c-baryon, the vertex is located 
closer to the horizon and consequently becomes less energetically favoured comparing to the 
iV-baryon. Similarly in the case of (N + k), even though not as obvious, adding k hanging 
and radial strings to the configuration of iV-baryon results in positive energy increase in the 
binding potential, making this configuration less favoured energetically. An (N + /c)-baryon 
naturally tends to decay into iV-baryon plus k free antiquark strings. A /c-baryon also has 
the tendency to fuse with (N — k) quarks to form an iV-baryon with lower energy. 

The situation of j-mesonance is somewhat similar. Even though j mesons are always 
energetically preferred over j-mesonance for all values of j, j-mesonance with higher value 
of j has stronger binding force than the lower ones as is shown in Fig. [31 From the energy 
viewpoint, j-mesonance will prefer to split into a number of j mesons. It is notable that the 
screening length of j-mesonance will approach the value of meson, L* meson) but it will never 
exceed L meson . 

For the case of (Af+fc)-baryon and j-mesonance, there exist the limits k — > oo and j — > oo. 
The first limit for (N + /c)-baryon leads to the zero-size configuration which saturates the 
zero-force condition. The second limit for j-mesonance leads to the mesonic limit where the 
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configuration is similar to the system of j mesons as we will see in the following. 

From Eqn. (1281) . since A(n) ~ (j/N)~ l , A(n) becomes negligible for large j/N. Therefore, 
we can neglect A(n) and obtain that Epi does not depend on j/N. Using asymptotic 
expansions, Eqn. (j2Tj) becomes 



£ - s I dy 



yn _ x n 



where 



y n - i 

e a) r i 1 - a c2,{n ~ 2) 

r( l_I) L2 /(n- 2 ) 

i?/ 2 r(i-i)r(i 



0(^ n ), (62) 



C(n) = 



2 - 



« r (|-^) 



■ 2 

Now, consider Eqn. (126jl . we find the screening length (half the distance between quarks 

at which the binding energy is zero) by setting E to t = 0. In the limit of j/N becoming very 
large, we can obtain L* from the condition 

S{L m ) = 0, (63) 

leading to 

r(i)r(i - -) i ( n_2 )/ 2 



2' 



C{n). (64) 

.«rWr( r -)J 

Again, the case n = 3 and n = 4 correspond to the Sakai-Sugimoto and the AdS-Schwarzschild 
gravity dual model respectively. This expression is exactly the same as the screening length 
of meson in the deconfined phase from Ref. [2offi It is no surprise since in the j — > oo limit, 
the hanging strings from the boundary exert force overwhelmingly, therefore the "weight" of 
the baryon vertex plus the tension of radial strings become negligible. Effectively, the end of 
hanging string at the vertex will feel zero force down and thus the slope u' c will be zero. As 
a result, the strings from the boundary will hang smoothly and appear similar to hanging 
strings in the case of the mesonic state. 

Even in the deconfined phase, we therefore perceive that in addition to free quarks, an- 
tiquarks, and gluons, there will also be mesons and multi-quark states. Due to the lower 
energy, there are more iV-baryons than {N + /c)-baryons and fc-baryons. The relative popu- 
lations can be estimated using the Boltzmann factor 

exp(-^), (65) 

determined by the corresponding binding energy E for each state. 

A more precise way of considering the deconfined phase is to use the grand canonical 
potential as the indicator for the stable phase. Following Bergman, Lifschytz, and Lip- 
pert [19], we consider three phases of the deconfined soup, a vacuum phase and a nuclear 



5 Our definition of the screening length is one- half of the definition in Ref. 
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phase with broken chiral symmetry, and a xS-QGP. For sufficiently high chemical potential 
and moderate temperature, the nuclear phase of the multiquark states is preferred over the 
vacuum and xS-QGP phase. Exotic nuclear states such as A;-baryon, (N + /c)-baryon, and 
j-mesonance are characterized by the number of radial strings n s hanging down from the 
D4-branes to the horizon. It is found that the multiquark states with n s > 0.5 are unstable 
thermodynamically. However, all of these exotic states with 0.5 > n s > are less stable 
than the normal iV-baryon with n s = 0. 

For each value of n s , there exists a triple point where the grand canonical potentials of 
the three phases are equivalent. Varying n s , this triple point will move along the phase 
transition line between the vacuum and the xS-QGP as is shown in Fig. [HJ The stable region 
of the nuclear phase shrinks as n s increases. As n s > 0.5, the nuclear phase becomes ther- 
modynamically unstable with respect to the density fluctuations for most of the parameter 
space. 

8 Conclusion 

The gravity dual picture of the deconfmed phase suggests that the binding energy or potential 
between quarks and antiquarks in this phase is nonzero due to the Coulombic piece of the 
interaction. Since the colorless condition is not required in the deconfmed phase, exotic 
configurations of the multiquark states are possible. We investigate three classes of these 
configurations, fc-baryon, (N + /c)-baryon, and j-mesonance. It is found that all of these 
configurations are less energetically favoured than the normal iV-baryon as well as being less 
stable thermodynamically. 

The dependence of the screening length on the parameters k, k, j is studied and the 
results are shown in Fig. HH6J The screening length of fc-baryon and j-mesonance are notably 
increasing with the values of k and j whereas the screening length of (N + fc)-baryon is a 
decreasing function of k. Interestingly, j-mesonance has saturated value of screening length 
equal to the screening length of meson as j ; — >• oo. 

The dependence on the quark mass of the binding potential at the leading order is derived 
and found to be ~ m 1_n (n = 3,4 for the Sakai-Sugimoto, AdS-Schwarzschild model). The 
linear quark-mass dependence of the rest energy that we naturally expect is included in the 
regulator and therefore not present in the binding potential. 

In order to consider phase diagram involving exotic nuclear phase, we consider the Sakai- 
Sugimoto model where the flavour branes D8 and D8 are introduced. The flavour D8-branes 
action is identified with the grand canonical potential of the relevant phase. The nuclear 
phase is considered in the limit when the D4-branes are pulled all the way up to the flavour 
branes. Exotic multiquark states with a number of strings stretched down to the horizon, i.e. 
n s > 0, become less stable than normal iV-baryon (n s = 0) since radial strings attached to 
the D4-branes pull the D4-D8 configuration closer to the horizon. Nevertheless, comparing 
to the vacuum and the xS-QGP phase, the nuclear phase of exotic multiquark states can be 
more stable in a region of phase diagram with high chemical potential and low temperature 
as is shown in Fig. [81 In this region, we expect to have a nuclear phase where iV-baryons, 
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fc-baryons, and (N + fc)-baryons coexist. For j-mesonance with n s = 1, our consideration 
of the grand canonical potential suggests that it is thermodynamically unstable to density 
fluctuations since ^ < 0. Generically, numerical studies reveal that exotic baryons with 
n s > 0.5 (namely /c-baryon with k/N < 0.5, (N + /c)-baryon with k/N > 0.5 and any j- 
mesonance) in the deconfined phase are thermodynamically unstable to density fluctuations. 

Acknowledgments 

We would like to thank Wen-Yu Wen, Ahpisit Ungkitchanukit and Kazuyuki Furuuchi for 
valuable comments. E.H. is supported by the Commission on Higher Education (CHE), 
Thailand under the program Strategic Scholarships for Frontier Research Network for the 
Ph.D. Program Thai Doctoral Degree for this research. P.B. and A.C. is supported in part 
by the Thailand Research Fund (TRF) and Commission on Higher Education (CHE) under 
grant MRG5180227 and MRG5180225 respectively. 



A Force condition at the D8-branes 



There are three forces acting on a D4 locating inside the D8-branes, one from the D8, 
another from the radial strings pulling down towards horizon and lastly the force from its 
own "weight" in the background. The equilibrium can be sustained only when these three 
forces are balanced. As is shown in Ref. [TJJ], variation of the total action with respect to u c 
and the constant of motion with respect to x±{u) lead to 



x A [u c ) 



L(u c ) 




k 2 (i + s 



where the Legendre transformed action is 



d 2 u~ 3 

7SJ 



5 



DS 



L(x' 4 (u), d) du, 
M / du u 4 ^/f(u)(x' 4 (u)) 2 + u- 3 



d 2 



ir' 



and the source term is given by 

*-> source. 



Md 



;UcVf( u c) +n s (u c 



(66) 
(67) 



(68) 
(69) 

(70) 



There are two contributions from the D-branes and strings as the sources for the baryon 
chemical potential. Additional strings increase the baryonic chemical potential of the exotic 



23 



multiquark states. Since the number of total charge on each D4 is N which is absorbed into 
TV, the number of radial strings stretched down to the horizon, n s , is thus given in unit of 
1/N. 
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